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Recent measurements of multijet production in neutral current deep inelastic ep scattering at HERA are presented. 
Emphasis is put on parton dynamics at low x. 
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5^; 1. INTRODUCTION 

o ; 

' Deep inelastic scattering (DIS) off protons has provided decisive information on the parton distribution functions 
(PDFs) of the proton. Inclusive measurements of the cross section for the reaction ep — > eX as a function of the 
virtuality of the exchanged boson (Q 2 ) and of the Bjorken scaling variable (x) have been used to determine the 
proton structure function Ff (x, Q 2 ). Perturbative QCD (pQCD) in the next-to-leading-order (NLO) approximation 
0^ ■ has been widely used to extract the proton PDFs from such measurements and to test the validity of pQCD. In the 
i— i standard approach (DGLAP [1]), the evolution equations sum up all leading double logarithms in lnQ 2 -lnl/x along 
with single logarithms in InQ 2 and are expected to be valid for x not too small. At low x, a better approximation 
is expected to be provided by the BFKL formalism [2] in which the evolution equations sum up all leading double 
logarithms along with single logarithms in In 1 jx. 

The DGLAP evolution equations have been tested extensively at HERA and were found to describe, in general, 
the data. In particular, the striking rise of the measured Ff (x, Q 2 ) at HERA with decreasing x can be accomodated 
in the DGLAP approach. Nevertheless, the inclusive character of Ff (x, Q 2 ) may obscure the underlying dynamics 
at low x, and more exclusive final states like forward 1 jets [3] need to be studied. BFKL evolution predicts a larger 
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fraction of small-x events containing high- Fx forward jets than predicted by DGLAP [3, 4]. Parton dynamics at 
low x is particularly relevant for the LHC given that most of the interesting hard processes involve partons with low 
fractional momenta. 
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oo 

9? 2. FORWARD JET PRODUCTION 

> : 

Forward jet production in neutral current (NC) DIS has been studied extensively at HERA. As an example, 
measurements of forward jet production withp t j ct > 3.5 GeV, polar angle 0j et between 7° and 20°, 0.5 < p 2 j Ct /Q 2 < 5 
and Xj c t = Fjot/Fp > 0.035, where E p is the proton-beam energy, in the kinematic region defined by 10~ 4 < x < 
4 • 10~ 3 and 5 < Q 2 < 85 GeV 2 are shown in Figure la and exhibit a strong rise towards low x [5]. Perturbative 
QCD apparently does not account for such a rise: the leading-order (LO) QCD calculation does not predict any 
rise while the one predicted by NLO [6] is still much too low (see Figure la). Some of the Feynman diagrams that 
are accounted for in the LO (0(ct s )) and NLO (0(ct 2 )) calculations are shown in Figure lb. The former has no 
additional gluon radiation and helps to understand why in the LO calculation there is hardly any phase space left 
for forward jet production. In contrast, the NLO calculation accounts for the radiation of one additional gluon. 
This explains why there is such a huge increase from LO to NLO: due to the opening of a new channel, namely 
gluon-exchange in the t channel. However, that means that the NLO calculation is effectively a "LO" calculation, 
since no corrections are included. The NLO calculation should thus have large theoretical uncertainties from higher 



1 The coordinate system used is a right-handed Cartesian system with the Z axis pointing in the proton beam direction, referred to as 
the "forward direction". 
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orders. A variation of the renormalisation scale around n 2 R = (p? dijets) °-° no ^ §i ve p i se to such large theoretical 
uncertainties (see Figure la). However, if Q 2 is instead chosen as the renormalisation scale, the resulting theoretical 
uncertainties are large, as pointed out in [5] and shown in Figure lc: measurements of forward jet production [7] with 
E^ > 5 GeV, pseudorapidity in the range 2 < ?7 jct < 4.3, 0.5 < (E^f/Q 2 < 2 and x iet > 0.036 in the kinematic 
region given by 4 • 10~ 4 < x < 5 • 10~ 3 and 20 < Q 2 < 100 GeV 2 are compared to NLO QCD calculations [6] with 
A*k = Q 2 ■ Large theoretical uncertainties which arise from higher orders in the pQCD calculations prevent a firm 
conclusion. Further progress can be made by making measurements for which genuine NLO calculations are available, 
i.e. one gluon radiation at LO and two additional radiated gluons at NLO. That is the case for three-jet production, 
which was already studied in [5] and has been investigated more thoroughly in [8]. The latter is discussed next. 
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Figure 1: Measurements of forward jet production in NC DIS as functions of x (a,c). Examples of Feynman diagrams (b). 



3. MULTIJET PRODUCTION AT LOW x 

Some of the Feynman diagrams accounted for in the pQCD calculations for three-jet production in NC DIS 
are shown in Figure 2a. A measurement of the differential cross section da/dx as a function of x for three-jet 
production [8] is presented in Figure 2b. The jets are reconstructed in the ~/*p frame and are required to fulfill the 
following conditions: the transverse momentum of each jet p t ^ > 4 GeV, the sum of the highest and next-to-highcst 
Pt jets above 9 GeV, the jet pseudorapidity in the laboratory frame to lie between -1 and 2.5 and at least one of the 
jets in the central region (—1 < T)^£ < 1.3). The kinematic region is defined by 10~ 4 < x < 10~ 2 , 5 < Q 2 < 80 GeV 2 
and 0.1 < y < 0.7, where y is the inelasticity variable. Perturbative QCD calculations are compared to the data in 
Figure 2b: the LO (0(a 2 )) calculation still falls short of the data, but the NLO (O(o^)) calculation [9] improves 
dramatically the description of the data at low x. 

The inclusion of the 0(a^) QCD corrections does not only improve the description of the measured rate but 
also that of the topology of the events. Measurements have been made of the distributions in the variables used 
to describe the topology of three-jet events in the three-jet centre-of-mass frame: the scaled energy of the jets 
X- = 2E'J(E[ +E' 2 + E' 3 ) (i = 1, 2; E[ > E' 2 > E' 3 ) and the two angles 9' and ip' (see Figure 2c). The measurements 
are shown in Figure 3 and compared to NLO calculations [9]. The inclusion of additional gluon radiation provides 
an improved description of the data, for example, in the cos 9' distribution: the NLO follows the data and exhibits 
peaks at cos 9' = — 1 and 1, whereas the LO calculation flattens out at cos 9' = — 1. 



2 



34 International Conference on High Energy Physics, Philadelphia, 2008 





0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

X', X' 2 cose' cosy' 

Figure 3: Measurements of the differential cross sections as functions of the variables used to describe the topology of three-jet 



events in the three-jet centre-of-mass frame. 



Further investigations of low- a; parton dynamics have been made by studying transverse-energy and angular cor- 
relations in dijet and trijet production in NC DIS [10]. Jets are reconstructed in the hadronic centre-of-mass (HCM) 
frame and required to fulfill the following conditions: -E^hom > ^ GeV, -E^'hcm > 5 GeV and —1 < rf^ ' ' < 2.5. 
The kinematic region is given by 10 -4 < x < 10~ 2 , 10 < Q 2 < 100 GeV 2 and 0.1 < y < 0.6. One of the most 
interesting angular correlations is provided by the variable \A<t>jlcM l> wn i cn is defined as the azimuthal separation 
of the two jets with largest -E^hcm- For dijet events, 0(a s ) kinematics constrain |A(/)^,^j 2 | to tt and 0(a 2 ) cal- 
culations provide the LO contribution; 0(ctg) calculations give the NLO correction. Measurements of the doubly 
differential cross section d 2 da / d\&(jy^}^ \dx for dijet production in different regions of x are presented in Figure 4. 
The 0(a 2 ) predictions increasingly deviate from the data as x decreases, whereas 0{a\) calculations [9] provide a 
good description of the data even at low x. 

In summary, parton dynamics at low x is vigorously pursued at HERA. Precise measurements of multijet production 
in NC DIS have been made down to x ~ 10 -4 in terms of jet rates, topologies and correlations. Comparison of 
perturbative QCD calculations with these measurements demonstrate the big impact of initial-state gluon radiation. 
Perturbative QCD at 0(a^) reproduces succesfully the measurements. However, the theoretical uncertainties are 
still significant and the precision of the data demands next-to-next-to-leading-order corrections to be included. 
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Figure 4: Measurements of the doubly differential cross section d da / dlA^Q-^ \dx for dijet production in different regions of x 
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